Gravity-Sensitive Quantum Dynamics in Cold Atoms 
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We subject a falling cloud of cold cesium atoms to periodic kicks from a sinusoidal potential created by a 
vertical standing wave of off-resonant laser light. By controllably accelerating the potential, we show quantum 
accelerator mode dynamics to be highly sensitive to the effective gravitational acceleration when this is close to 
specific, resonant values. This quantum sensitivity to a control parameter is reminiscent of that associated with 
classical chaos, and promises techniques for precision measurement. 
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The identification and observation of signatures of chaos in 
quantum dynamics is the goal of considerable current effort. 
Much of this work centers on the theoretical definition and 
characterization of energy spectra [1], or such quantities as 
the Loschmidt echo J3l and fidelity [3], which essentially de- 
velop the idea that sensitivity of a wavefunction's evolution to 
small variations in a system's Hamiltonian be used as a defini- 
tion of quantum instability 111 H yl tHI ■ Such quantities could 
be observed experimentally but require some interpretation to 
highlight the way in which their nature betokens stability or 
chaos. An attractive alternative would be the observation of 
different motional regimes. This is more in sympathy with 
the techniques and philosophy used to identify classical chaos, 
and is the approach used here. 

In certain systems the decay of the overlap of two ini- 
tially identical wavefunctions evolving under slightly differing 
Hamiltonians can be expressed in the long time limit as the 
sum of two decay predictions, governed by Fermi's Golden 
Rule and the classical Lyapunov exponent |2j. The decay 
rate serves as a quantum signature of instability, which can 
be compared with that of the corresponding classical system. 
Such sensitivity can be probed by interferometric techniques 
d^l, |^] . In the quantum-mechanical system presented here, the 
classical limit of which is chaotic, extreme sensitivity of the 
qualitative nature of the motional dynamics to a control pa- 
rameter is directly observable. It is manifested by the effect on 
quantum accelerator mode (QAM) dynamics 0000 Holl 
of small variations in the effective value of gravity in the 6- 
kicked accelerator an extension of the paradigmatic 6- 
kicked roto r 1 11 ll The QAM observed in this atom optical 
realization are characterized by a momentum 

transfer, linear with kick number, to a substantial fraction (up 
to ~ 20%) of the initial cloud of atoms. This is due to a reso- 
nant rephasing effect, dependent on the time-interval between 
kicks, for certain initial wavefunctions Ii9l 1 1 211 - The sensitiv- 
ity in the dynamics we observe also promises the capability 
of precisely calibrating a relationship between the local grav- 
itational acceleration and /i/m, where m is the atomic mass, 
and we describe how our observations constitute a prelimi- 
nary feasibility-demonstration of such a measurement. 

The Hamiltonian of the ^-kicked accelerator, realized using 



a magneto-optic trap (MOT) of laser-cooled atoms that are 
then released and subjected to pulses from a standing wave of 
off-resonant light, is 

H= E-+mgz- h(p d [\ + cos(Gz)] V S(t - nT), (1) 
2m 

« = -oo 

where z is the position, p the momentum, m the particle mass, 
t the time, T the pulse period, ntpd quantifies the strength of the 
kicking potential, G = 27r//i spa t, and /l spat is the spatial period 
of the standing wave applied to the atoms. The quantity g is 
normally the gravitational acceleration. However, by 'acceler- 
ating' the standing wave, it is possible to effectively modify g. 
We have previously used this technique to counteract gravity 
and regain kicked rotor dynamics I7l ll3ll . 

In an innovative analysis by Fishman, Guarneri, and Re- 
buzzini (FGR) d, the fact that QAM are observed only 
when T approaches fT\/2 = f2mn/HG 2 , where I e Z + and 
T\j2 is the half-Talbot time |9], is exploited to yield a dra- 
matically simplified picture of QAM dynamics. In a frame 
accelerating with g, the linear potential is removed to leave a 
spatially periodic Hamiltonian. The quasimomentum f3 is then 
conserved, i.e., if a momentum state \p) - \(k + f3)hG), where 
k e Z and f3 e [0, 1), 'ladders' of momentum states of dif- 
ferent j3 evolve independently. The resulting «-dependent and 
yS-specific kick-to-kick time evolution operator is 



F n (J3) = exp(-/{p + sgn(e)[7rf + k/3 - y(n - l/2)]} 2 /2e) 
x exp(/fccos;£-/|e|), 



(2) 



where k = \e\<pd, k = 2nT/Ti/2, and y = gGT 2 . We have in- 
troduced a smallness parameter, e = 2n(T/Ti/2 - (), to quan- 
tify the closeness of T to {T1/2 and the dynamical variables 
are now an angle x - Gz and a discrete conjugate momen- 
tum p = p\e\/hG, such that \x,p] = i\e\. If one constructs a 
kick-to-kick Heisenberg map corresponding to Eq. for the 
dynamical variables, then in the limit e — > 0, the commutator 
vanishes along with the uncertainty principle, and the opera- 
tors can be replaced by their mean values. Thus 



Pn+i = Pn - ksin(x„) - sgn(e)y, 
Xn+\ = Xn + sgn(e)p„+i, 



(3a) 
(3b) 
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FIG. 1: Phase space plots produced by Eq. 0, when Iny/k 1 = 
1/1 => y = (2tt + e) 2 /2?r, and £ = |e|0.8w for (a) e = -0.88, 
(b) 6 = -0.02, (c) 6 = 0.03, (d) 6 = 0.6. This corresponds to 
T = 57.4 fis, 66.5 //s, 67 ^s, 73 //s. For (a), (b) the island corresponds 
to a (p, j) = (1, 1) QAM, and for (c), (d), to a (p, j) = (1, -1) QAM. 



where 6„ = 0„) and p„ = (p n ) + sgn[7r£ + k/3 - y(n - 1/2)]. 
Quantum accelerator modes are one-to-one related to stable 
periodic orbits of this map 0, [H O . It is very important to 
note that e — » coincides with h — > only if £ = 0. Oth- 
erwise, as in the experiments here, the classical-particle-like 
behavior of QAM is due to a quantum resonance effect. 

The stable periodic orbits yielded by Eq. (0 (and hence 
QAM) are classified by their order p and jumping index j (the 
number of momentum units, in terms of the size of the phase- 
space cell, traversed after p iterations). The sign of j is deter- 
mined by whether this is in the positive or negative momen- 
tum direction. A necessary condition 11211 for the existence of 
a periodic orbit is |j/p + sgn(e)y/2^| < k/2n, which can be 
rewritten (for small e) as 



\2n 



(4) 



Both 4>d an d y/k = gm /fr G are independent of T, and 
therefore of e. Equation is convenient when T is varied 
from just below to just above £T\/2, i.e., scanning e from neg- 
ative to positive, as in the experiments described here. As 
e — > 0, the QAM that occur must be characterized by j and 
p such that j/p — > -sgn(e)27rf 2 y/£ 2 . In general Infiy/k 2 is 
an irrational value, and one usually observes a succession of 
increasingly high-order QAM as T — » (T\j2 [;83. If we tune 
g so that 2n( 2 ylk 2 — r/s, where r and s are integers, then 
j/p + sgn(e)27rf 2 y/£ 2 = for j/p = -sgn(e)r/s. Once the (p, j) 
QAM satisfying this condition appears, shifting T closer to 
{Ti/2 does not result in higher-order QAM. 

In a frame accelerating with g, the momentum after 



kicks, for an initial condition near a (p, j) stable periodic orbit 
Jl2ll . in 'grating recoils' fiG is 



2tt 

q N - q + N— 



- +sgn(e) — 

P 27T 



(5) 



where go is the initial momentum. For jV a multiple of j, this 
result is exact for e-classical initial conditions located on (p, j) 
periodic orbits. We now consider the momentum of orbits 
specified by j/p = r/s (for e < 0) and j/p = —r/s (for e > 0) 
as a single function of and e, in the case where 2n( 2 y/k 2 
approaches rational values. Letting 2n( 2 ylk 2 — r/s + w€ 2 
0, we find: 



r 12 e 
qN ^ qo + N -\- + — 



<2n( 2 e \ 

Nw I +2(+—\. (6) 

e 27T/ 



Scanning through e from negative to positive values, one does 
not generally observe two QAM of the same p and magnitude 
of j (with positive sign for negative e, and negative sign for 
positive e) I8l ll5ll . However, in the gravity -resonant cases we 
consider, when 2n( 2 y/k 2 is close to r/s, we always observe 
an r/s and then a -r/s QAM as we scan e in this way. This 
is shown in Fig.^ where we plot Poincare sections produced 
by Eq. Q for y = k 2 /2n = (2n + e) 2 /2n (i.e., r/s = 1) and 
k — \e\<pd — |e|0.87r (the approximate experimental mean value 
0]). The islands around the (p, j) = (l,-sgn(e)l) periodic 
orbits remain large over a wide range of e and, in dramatic 
contrast to Ref. |8], no higher-order island structures appear 
as e — > 0. In Fig. [2 a) the corresponding QAM are similarly 
robust and uninterrupted by higher-order QAM as T — > 7\/2. 

From Eq. we thus see that for a given N, q is a linear 
function of e whenever w = 0. If w £ this changes to a 
hyperbolic function of e, where the arms of the hyperbolae 
point in opposite directions for oppositely signed w. Devi- 
ation from straight line behavior in a QAM accelerated to a 
given momentum will be greater for a gravity-resonant mode 
corresponding to a smaller value of j/p = rj s. This is because 
the acceleration of the mode is oc j/ p, but the deviation is oc N. 
We consider only QAM where j = r = 1, so high-order modes 
exhibit, for a given momentum transfer, greater sensitivity to 
variations in g than low-order modes. 

In our realization of the quantum (5-kicked accelerator, 
~ 10 7 cesium atoms are trapped and cooled in a MOT to a 
temperature of 5fiK, yielding a Gaussian momentum distribu- 
tion with FWHM 6hG. The atoms are then released and ex- 
posed to a sequence of equally spaced pulses from a standing 
wave of higher intensity light 15 GHz red-detuned from the 
6 2 Si/ 2 -> 6 2 Pi/ 2 , (F = 4 -> F' = 3) Dl transition. Hence 
the spatial period of the standing wave is A spat = 447 nm, and 
7*1/2 = 66.7 [is. The peak intensity in the standing wave is 
- 5 x 10 4 mW/cm 2 , and the pulse duration is t p = 500 ns. 
This is sufficiently short that the atoms are in the Raman- 
Nath regime and hence each pulse is a good approximation 
to a (5-function kick. The potential depth is quantified by 
<pd = ^ 2 fp/8(5i, where Q. is the Rabi frequency and 5l the 
detuning from the Dl transition. During the pulse sequence, 
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FIG. 2: (color online). Color density plots of experimental momentum distributions for different effective gravity g corresponding to (a) 
r/s = 1/1 (after 15 kicks), (b) r/s = 1/2 (30 kicks), (c) r/s = 1/3 (45 kicks), and (d) r/s = 1/4 (60 kicks), as T is varied in the vicinity of 
the half Talbot time Tip = 66.7 /js, from 60.5 fis to 74.5 //s in steps of 0.128 fjs. In each case the QAM corresponds to j/p = r/s; subplot (i) 

corresponds to w - -8.5 x 10~ 4 (deviation from resonant g is 8.6 x 10~ 2 ms~ 2 ), subplot (ii) to w =s 0, and subplot (iii) to w - 8.5 x 10~ 4 

(deviation from resonant g is ~ 8.6 x 10 2 ms~ 2 ). Overlaid lines, labeled (p, j), indicate QAM momenta predicted by Eq. |6j. Population 
arbitrarily normalized to maximum value = 1, and momentum denned in a frame falling with g. Note the significantly greater population at 
high momentum (up to 50UG) near T1/2 in (d.i) and (d.iii), compared to (a.i) and (a.iii). 



a voltage-controlled crystal phase modulator is used to stro- 
boscopically accelerate the standing wave profile. The atoms 
therefore effectively experience a non-standard, and control- 
lable, value of gravity. After the pulsing sequence, the atoms 
fall through a sheet of laser light resonant with the 6 2 S 1/2 — > 
6 2 P 3/2 , (F = 4 -> F" = 5) D2 transition, 0.5 m below the 
MOT. By monitoring the absorption, the atoms' momentum 
distribution is then measured by a time-of-flight method, with 
resolution tiG. For further details see Refs. ||7|,|9||. 

In Fig. 13 we show momentum distributions for experiments 
in which the value of T was scanned around 7*1/2 (( = 1) 
from 60.5//S to 74.5//S, with 27ry/^ 2 varied in the vicinity of 



r/s equal to (a) 1/1, (b) 1/2, (c) 1/3, and (d) 1/4. To main- 
tain the ideal (w = 0) total momentum transfer, 15, 30, 45 
and 60 kicks were applied, respectively, fixing Nr/s. For each 
of Figs. |3 a), |3Jb), 0c), and|5Jd) the data displayed are: in 
subplot (ii), from experiments in which 27ry/£ 2 = r/s is ful- 
filled as exactly as feasible, yielding linear variation of the 
QAM momentum with T; and in subplots (i) and (iii), for 
equal positive and negative deviations, respectively, from this 
near-ideality, yielding hyperbolic variation of the QAM mo- 
mentum. Typically ~ 10-20 % of the atoms are accelerated 
away from the cloud centered at p = 0. 

In each subplot (ii) of Fig. |2j the QAM momentum pre- 
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dieted by Eq. 0, shown as an overlaid line, is identical. The 
expected linear dependence on T appears to be well confirmed 
by the data, although the separation of the QAM from the 
main, non-accelerated cloud, centered at zero momentum, is 
clearer for smaller s — p (there is less momentum diffusion 
of the main cloud due to the smaller number of kicks). The 
effect of imperfectly resonant gravity, shown in subplots (i) 
and (iii) in Fig. 13 is much more dramatic for larger s = p, for 
which more kicks are applied. In Fig. EJa), subplots (i) and 
(iii) are barely distinguishable from subplot (ii), whereas in 
Fig.|3d), the momentum distributions in subplots (i) and (iii) 
are highly asymmetric compared with subplot (ii), with, close 
to T\/2, noticeable population at up to 50hG. The asymmetry 
inverts as one changes from below [subplot (i)] to above [sub- 
plot (iii)] the resonant value of gravity. We therefore observe 
a clear qualitative change in the QAM dynamics, highly sen- 
sitive to a control parameter. The displayed predictions of Eq. 
(|6} show that deviations from linear behavior only occur when 
very close to T\/2 in Figs.|2ja.i) and|2ja.iii), but are much more 
significant in Figs.|2jd.i) and|2[d.iii). This is due to the larger 
number of kicks necessary for large s = p to achieve the same 
QAM momentum. 

The procedure of determining the 'standing wave acceler- 
ation' at which straight-line behavior of a given (p, j) QAM 
momentum is observed as a function of T could, in principle, 
be used as a sensitive atom-optical means of relating h/m lllfill 
to the local gravitational acceleration lll7ll . This is because 
Iny/k 1 — r/s can be rephrased as g — (h/m) 2 (r/s)/ t and 
would be determined by noting when the total acceleration 
(sinusoidal potential plus gravitational) causes these equali- 
ties to be fulfilled for a known rj s, and then subtracting the 
imposed acceleration of the potential. In our setup, where the 
sinusoidal potential is 'accelerated' by using a crystal phase 
modulator to phase-shift the retroreflected laser beam |7l|9|], 
the value of the phase shift due to a particular applied volt- 
age is difficult to calibrate more precisely than ~ 1 %. This 
accordingly limits our measured precision of the relationship 
between the local gravitational acceleration and h/m to ~ 1 %. 
Accurate prediction of the QAM momenta for imperfectly res- 
onant values of the effective gravity, as displayed in Fig.|2] is 
also hampered. This could be improved by a configuration in 
which a moving sinusoidal potential is formed by two coun- 
terpropagating beams with a controllable frequency difference 
111 (811 . where calibration of the phase shift to between lppm and 
lppb is possible. Calibration of /i spat to less than lppb is also 
feasible fl7ll . allowing for the possible sensitive determina- 
tion of either the local gravitational acceleration [17] or h/m 
lUrill . depending on which is known more precisely at the out- 
set. The feasibility of any such scheme will ultimately depend 
on how precisely the atomic ensemble's dynamics permit the 
determination of the acceleration of the sinusoidal potential 
for which the resonant, linear with T, behavior of the QAM 
occurs. Ascertaining this will require substantial theoretical 
and experimental investigation. 

In conclusion, we have observed qualitative changes in the 



motional quantum dynamics of cold cesium atoms, which are 
highly sensitive to the precise value of an externally adjustable 
parameter, the effective gravity. This is distinct from concep- 
tually related proposals that consider slightly differing Hamil- 
tonians to study the Loschmidt echo or fidelity, and demon- 
strates an attractive link to the concepts of highly sensitive dy- 
namics in classically chaotic systems. Furthermore, we have 
described a feasible experimental scheme taking advantage of 
this sensitivity to determine a relationship between the local 
gravitational acceleration and h/m. 
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